IN THE CLAIMS : 

1. (Original) An image reconstruction method comprising: (a) obtaining the 
measurement data corresponding to the image; (b) processing said measurement data using 
heuristic non-linear global optimization , specifically the Method of Simulated Annealing in 
order to obtain an image; and (c) displaying said processed image on a display device. 

2. (Original) A method according to claim 1, wherein the measurement data are 
electrical parameters on the perimeter of a region such as the inside of a pipeline, of the oil 
well or of the a tank, for the obtained image of the spatial distribution of the electrical 
permittivity (or dielectric constant) or of the electrical conductivity within said region, and 
the image reflects the spatial distribution of the materials or substances that occupy the 
region, such as gases or/and liquids. 

3. (Original) A method according to claim 2, wherein said parameters are the 
electrical capacitance values measured between the electrodes of a sensor, which is formed 
by a plurality of said electrodes placed around the perimeter of the region (pipeline, well, 
tank). 

4. (Currently amended) A method according to claims 2 and 3 claim U wherein said 
the image obtained is an image of the spatial distribution of the electrical permittivity (or 
dielectric constant) within said interest region (pipeline, well, tank)., which reflects the 
spatial distribution of the materials or substances, such as gases or/and liquids, that occupy 
this region (pipeline, well, tank). 



5. (Currently amended) A method according to claimo 2 to ^ claim 2. wherein said 
the image is formed by a finite number of sub-regions or pixels in the viewing, which number 
depends on the desired resolution. 



6. (Original) A method according to claim 3, wherein said sensor is formed by a 
pipe made of an electrically insulating material, on whose outer wall an array of rectangular 
metallic electrodes is placed. 

7. (Currently amended) A method according to claims 3 and 6 claim 3 , wherein the 
sensor contains a multiphase or multicomponent flow, and images that show the distribution 
of the phases or components, such as gases or/and liquids, are obtained in accordance to what 
is said in the background of the invention. 

8. (Original) A method according to claim 1, wherein the method of simulated 
annealing is used as the optimization method. 

9. (Original) A method according to claim 8, wherein a cost function associated to 
the energy of the system is iteratively minimized with respect to e (the vector of the 
permittivity values in each pixel of the image), said cost function being of the form: 



where cT^""^ are the m measured mutual capacitances and cf''''(e^) are the ones calculated by 

solving the forward problem for a given permittivity distribution Zk - 
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10. (Original) A method according to claim 9, wherein the Metropolis criterion is 
used to do the minimization. 

11. (Original) A method according to claim 10, wherein in the minimization process 
it is used, as the initial guess for the permittivity distribution, 8, both a homogeneous 
distribution and the distribution that results from applying the linear back-projection (LBP) 
method to the measurement data. 

12. (Original) A method according to claim 9, wherein the computation of the 
calculated capacitances cf'''''(Ej^) , known as the forward problem, is carried out by means of 
the finite- volume method. 

13. (Original) A method according to claim 12, wherein for the solution of the 
system of equations that results from solving the forward problem cf''^''(ef^), iterative 
methods that rapidly converge are used, by employing as initial estimation of the electrostatic 
potential the result of the forward problem solution obtained in the previous iteration of the 
inverse problem. 

14. (Original) A method according to claim 1, wherein the method of genetic 
algorithms is used as the optimization method. 



15. (Original) A method according to claim 14, wherein, starting from an initial 
population of Q permittivity models Bk{o) {k= I, ,„,Q\ evolutionary mechanisms such as 
selection, crossover and mutation are applied in order to obtain new populations. 



16. (Original) A method according to claim 15, wherein individuals 8* are 
characterized by having a small cost (or misfit) function, which is given by 



= ^lik) = jrz : (/ - 1, . . , m) 



where c"^" are the m measured mutual capacitances and cf^^'^iBf^) are the calculated ones by 
solving the forward problem for a given model or permittivity distribution . 

17. (Currently amended) A method according to claims 15 and 16 claim 15. wherein 
the accumulated probability of selection for a particular model is given by 

^ \ max avr y 

where Emax, and Eavr are maximum and average cost functions of the generation, respectively, 
and Q is the number of individuals in the population. 

18. (Original) A method according to claim 17, wherein a biased roulette procedure 
can be used to decide which models are selected on each iteration of the method. 



19. (Original) A method according to claim 18, wherein crossover and mutation of 
models are carried out randomly according to the probabilities of crossover and mutation, Pc 
yPm^ 



20. (Original) A method according to claim 19, wherein the probability of mutation 
Pm is determined using the average variation coefficient y, given by 



1 p 

p 1^ 



r \ 



where p is the number of parameters, £. is the average of the /-th parameter, and o; is the 
standard deviation. 

21. (Currently amended) A method according to claims 19 and 20 claim 19 , wherein 
Pni is defined as a function of /, that is: 



/L = 



Pini para ;^>0.1 

0.1 para 0.02 <:k< 0.1 

0.2 para y < 0.02 



where Pi„i is the initial probability of mutation. 



22. (Currently amended) A method according to claims 15 and 16 claim 15. wherein 
the calculation of cf'''''(Si^), known as the forward problem, is performed by means of the 
finite- volume method. 
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